Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1436 1479 (N30W33 L= 22)   240     33        1 (1415/1460-1459/1500) -56.3 W31-W35 N28-N31
      2 1584 1364 (N22W42 L= 30)   256     58        0 (1571/1331-1599/1399)  76.0 W40-W43 N19-N23
      3  310 1070 (N09E50 L=299)   304    100        0 ( 300/1025- 339/1110)  59.0 E50-E47 N06-N11
      4  154 1003 (N05E68 L=280)  3413    211        2 ( 122/ 913- 202/1105)  84.8 E74-E61 N00-N11
      5  443  981 (N04E38 L=310)   792    160        2 ( 411/ 909- 473/1068)  77.5 E41-E35 N00-N09
      6  933  925 (S01E06 L=342)   388    145        0 ( 887/ 874- 956/ 991)  45.0 E09-E04 S03-N02
      7  279  758 (S10E56 L=293)  1759    285        4 ( 160/ 665- 391/ 819) -16.3 E70-E46 S15-S06
      8 1097  716 (S15W03 L=351)   452    139        2 (1039/ 664-1139/ 747) -31.0 W00-W05 S18-S12
      9  354  592 (S21E53 L=295)   728    147        1 ( 291/ 557- 424/ 629) -18.4 E59-E46 S23-S18
  

    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