Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  662 1617 (N51E08 L=265)   401     92        0 ( 645/1573- 678/1662)  82.9 E09-E06 N47-N55
      2 1335 1602 (N29W38 L=311)   409     48        1 (1303/1580-1369/1618) -21.8 W36-W41 N26-N30
      3  419 1171 (N28E36 L=237)   484     81        1 ( 377/1134- 453/1209) -45.0 E38-E32 N24-N30
      4  500 1082 (N20E31 L=242)   413     88        0 ( 459/1053- 547/1102)   8.1 E33-E26 N19-N21
      5  212  927 (N17E59 L=214)   391     67        1 ( 193/ 898- 235/ 960)  36.9 E61-E55 N15-N18
      6 1868  786 (S34W58 L=331)  1557    124        4 (1838/ 727-1896/ 854)  76.0 W53-W63 S37-S29
      7 1557  649 (S32W24 L=297)   463    123        0 (1510/ 611-1599/ 702)  20.6 W19-W28 S33-S29
      8 1707  645 (S37W38 L=311)   430     51        1 (1682/ 622-1736/ 674)  45.0 W35-W41 S39-S35
      9 1208  549 (S28E02 L=271)   343    136        0 (1154/ 504-1261/ 590)  14.0 E06-W02 S29-S25
     10 1343  229 (S59E04 L=269)   557     61        2 (1312/ 212-1374/ 252) -50.2 E07-E01 S62-S56
  

    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