Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  513 1589 (N31E48 L=340)  1169    192        1 ( 450/1517- 569/1675) -26.6 E61-E41 N26-N36
      2  344 1454 (N21E57 L=331)   780    201        0 ( 278/1404- 397/1538)  68.2 E63-E54 N16-N26
      3  303 1500 (N23E65 L=323)   583     38        1 ( 281/1472- 317/1524) -45.0 E69-E61 N21-N24
      4 1055 1257 (N16E00 L= 28)   238     61        1 (1033/1231-1088/1289) -56.3 E01-W01 N14-N17
      5 1033 1057 (N03W00 L= 28)   489    185        0 ( 962/1016-1113/1089)  14.0 E03-W05 N00-N05
      6 1427  741 (S12W29 L= 57)  1617    473        1 (1260/ 633-1568/ 913)  51.3 W19-W36 S20-N00
      7  224  687 (S28E63 L=325)   656     90        0 ( 197/ 652- 263/ 718)  23.2 E70-E56 S31-S25
  

    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