Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  367 1594 (N41E57 L=333)   501     93        0 ( 341/1575- 397/1608)  26.6 E59-E54 N39-N42
      2 1681 1547 (N34W56 L= 86)   844    111        2 (1642/1506-1708/1579) -68.2 W52-W58 N31-N36
      3  536 1480 (N34E36 L=354)   954    232        4 ( 458/1429- 622/1523) -11.3 E43-E28 N30-N36
      4  375 1284 (N20E46 L=344)  2771    683        5 ( 233/1081- 476/1436) -36.9 E64-E36 N07-N29
      5  947 1386 (N27E04 L= 26)   372     96        1 ( 886/1372- 998/1402)   0.0 E08-W00 N26-N27
      6 1167 1179 (N13W09 L= 39)   264     83        0 (1124/1143-1209/1207) -38.7 W06-W11 N10-N15
      7 1563  990 (N00W35 L= 65)   364     80        1 (1520/ 971-1605/1013)   0.0 W31-W38 N00-N01
      8  496  913 (S02E35 L=355)   309     69        0 ( 465/ 898- 531/ 928)  11.3 E37-E32 S03-S01
      9  788  832 (S07E15 L= 15)   396    113        0 ( 726/ 799- 839/ 863) -29.7 E19-E12 S09-S05
     10  376  677 (S18E47 L=343)  1073    311        3 ( 314/ 580- 429/ 780)  47.1 E55-E41 S24-S10
     11 1436  580 (S26W28 L= 58)   585    140        1 (1375/ 566-1495/ 597)   0.0 W23-W32 S26-S24
     12  603  532 (S27E32 L=358)   996    250        2 ( 535/ 491- 706/ 594)  26.6 E38-E22 S30-S22
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y