Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  524 1553 (N34E42 L= 10)  1355    326        3 ( 429/1466- 592/1621) -22.8 E52-E33 N28-N38
      2 1068 1321 (N21W02 L= 54)  1239    298        0 ( 925/1235-1200/1427) -26.6 E08-W10 N16-N27
      3  221 1354 (N18E65 L=347)  1074    171        2 ( 198/1304- 245/1407) -80.5 E68-E60 N15-N21
      4  301 1330 (N18E55 L=357)   421    153        0 ( 260/1292- 355/1361)  26.6 E58-E50 N15-N20
      5  581 1250 (N15E30 L= 22)   588    172        1 ( 563/1172- 604/1335) -78.7 E32-E27 N09-N19
      6  733 1272 (N17E20 L= 32)   346    109        1 ( 701/1222- 769/1327) -60.3 E22-E17 N13-N20
      7  373 1127 (N06E45 L=  7)  1444    294        4 ( 330/1026- 441/1260) -52.4 E50-E38 N00-N13
      8  304  993 (S03E50 L=  2)   436    158        0 ( 254/ 978- 365/1012)  -5.2 E55-E44 S03-S01
      9 1091  266 (S51W11 L= 63)   617    131        1 (1033/ 229-1137/ 294) -28.6 W05-W16 S54-S48
  

    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