Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1218 1322 (N21W18 L=195)   474    192        0 (1138/1301-1304/1343) -21.8 W13-W23 N18-N23
      2 1408  928 (S06W21 L=199)   208     79        1 (1393/ 898-1426/ 970) -90.0 W20-W22 S08-S03
      3 1830  873 (S20W56 L=233)   791    141        0 (1802/ 806-1862/ 934)  45.0 W52-W60 S24-S17
      4 1686  892 (S15W41 L=219)   343     80        0 (1650/ 858-1710/ 929)  26.6 W37-W43 S16-S12
      5 1782  664 (S33W52 L=230)   835    141        0 (1754/ 633-1820/ 718) -54.5 W48-W58 S35-S28
      6 1545  583 (S32W28 L=205)  1609    280        2 (1484/ 482-1626/ 707)  45.0 W22-W34 S38-S25
      7 1051  648 (S16E05 L=172)   678    266        1 ( 983/ 605-1136/ 701) -20.6 E09-W00 S18-S12
      8 1013  382 (S34E15 L=163)  3059    453        5 ( 884/ 226-1154/ 543)  32.8 E32-E01 S45-S24
      9  542  518 (S17E43 L=135)  1469    345        3 ( 384/ 487- 695/ 542)  -9.1 E56-E30 S19-S13
     10  865  237 (S44E34 L=143)   380     44        1 ( 845/ 210- 885/ 266)  56.3 E36-E32 S46-S40
  

    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