Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1181 1615 (N33E01 L=341)   414     80        2 (1134/1598-1223/1630)  31.0 E04-W01 N31-N34
      2 1052 1495 (N22E07 L=335)   300     70        0 (1022/1468-1087/1520) -11.3 E09-E04 N20-N23
      3  994 1423 (N16E09 L=333)   223     64        0 ( 965/1407-1033/1448) -14.0 E10-E06 N15-N17
      4 1632 1237 (N18W35 L= 17)   243     82        0 (1606/1201-1652/1276) -45.0 W32-W36 N15-N19
      5 1809 1098 (N14W52 L= 34)   421     78        0 (1793/1058-1827/1140) -45.0 W50-W54 N11-N15
      6  289  962 (S20E50 L=292)   945    214        0 ( 231/ 893- 343/1040)  42.0 E56-E45 S24-S15
      7 1178  675 (S24W16 L=358)   306     75        1 (1141/ 646-1215/ 703)  68.2 W14-W18 S26-S21
      8 1634  628 (S17W49 L= 32)  1183    271        5 (1511/ 581-1738/ 671)  27.9 W40-W57 S21-S12
      9  960  582 (S34W04 L=347)   264    108        0 ( 904/ 568-1010/ 599)  23.2 W00-W07 S35-S32
     10  940  275 (S57W14 L=356)   630    101        1 ( 883/ 251- 984/ 306)   4.8 W08-W20 S58-S54
  

    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