Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1335 1571 (N32W36 L= 29)   505     99        1 (1284/1550-1385/1587) -15.9 W31-W39 N30-N33
      2  693 1454 (N40E14 L=340)   607    105        0 ( 639/1407- 747/1492)  16.7 E18-E08 N37-N41
      3 1044 1247 (N20W06 L=360)   644    200        3 ( 957/1226-1158/1272) -26.6 W00-W12 N16-N22
      4 1335 1195 (N10W22 L= 16)   649    165        1 (1244/1160-1404/1212) -24.0 W17-W26 N07-N13
      5  528  645 (S07E39 L=314)   286     29        0 ( 510/ 623- 550/ 670) -71.6 E40-E37 S07-S04
      6 1017  577 (S20E10 L=344)   429    146        1 ( 983/ 510-1053/ 634)  45.0 E13-E06 S23-S16
      7  507  535 (S13E45 L=309)   488     84        0 ( 486/ 496- 541/ 580) -90.0 E46-E43 S16-S10
      8  860  312 (S35E30 L=324)   467     59        2 ( 834/ 286- 883/ 336)  31.0 E32-E27 S37-S33
      9 1268  250 (S53E02 L=351)   369     54        0 (1240/ 238-1299/ 261) -45.0 E04-W00 S54-S50
  

    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