Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  711 1539 (N39E10 L=109)   372    112        1 ( 661/1520- 776/1556) -12.5 E14-E05 N37-N41
      2 1194 1415 (N21W18 L=138)   539    163        3 (1141/1386-1254/1439) -14.0 W14-W22 N19-N22
      3  444 1385 (N35E33 L= 86)   586    164        1 ( 388/1348- 481/1434) -41.6 E38-E29 N31-N39
      4  709 1155 (N15E16 L=103)   284     55        1 ( 686/1124- 737/1185) -90.0 E17-E14 N13-N17
      5 1321  969 (S07W16 L=135)   879    296        1 (1179/ 910-1448/1016) -20.6 W07-W23 S13-S03
      6 1067  791 (S13E02 L=117)  1181    335        2 (1010/ 627-1174/ 923)  49.1 E07-W06 S21-S06
      7 1088  596 (S25E05 L=114)   291     85        1 (1042/ 579-1133/ 623) -31.0 E07-E02 S26-S22
      8 1066  386 (S39E12 L=107)   240     64        1 (1027/ 372-1099/ 397)   0.0 E15-E09 S40-S37
  

    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