Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1512 1747 (N31W59 L=201)  2178    429        7 (1419/1696-1610/1784) -11.3 W49-W69 N27-N32
      2 1305 1605 (N29W36 L=178)   279     67        0 (1276/1589-1333/1625)   0.0 W33-W38 N28-N30
      3 1059 1060 (N06W03 L=145)   299    108        0 (1025/1003-1110/1111)  33.7 W00-W06 N03-N07
      4 1223 1038 (N00W11 L=153)   401    117        0 (1168/ 986-1274/1072) -63.4 W08-W13 S03-N03
      5 1859  917 (S25W54 L=196)   533    103        0 (1830/ 871-1881/ 963)  26.6 W50-W56 S27-S22
      6 1347  920 (S10W15 L=158)   207     87        0 (1307/ 905-1398/ 929) -18.4 W12-W18 S10-S08
      7 1606  720 (S29W29 L=171)   689    170        0 (1551/ 635-1660/ 805)  38.7 W23-W33 S32-S24
      8 1244  603 (S25W02 L=144)   953    310        1 (1150/ 483-1360/ 728)  30.5 E06-W11 S30-S19
      9  503  357 (S22E56 L= 86)  2729    525        4 ( 456/ 227- 534/ 503)  51.3 E66-E47 S32-S12
     10  587  442 (S18E44 L= 98)   304     38        0 ( 567/ 417- 604/ 461) -90.0 E45-E42 S20-S16
     11  885  383 (S30E28 L=114)   414    123        2 ( 867/ 333- 903/ 434)  60.3 E30-E25 S33-S26
     12  663  355 (S26E44 L= 99)   763    179        2 ( 635/ 289- 682/ 417)  48.4 E47-E39 S30-S21
  

    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