Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1218 1371 (N17W02 L= 49)   204     62        0 (1192/1332-1238/1403) -36.9 W00-W04 N15-N18
      2  321 1178 (S13E49 L=357)  1346    191        3 ( 272/1080- 386/1238)  63.4 E54-E44 S19-S08
      3  683 1025 (S15E20 L= 26)   544    199        0 ( 632/ 937- 735/1113) -38.7 E24-E14 S18-S10
      4 1367  958 (S02W21 L= 67)   989    233        1 (1270/ 905-1473/1035)   0.0 W13-W27 S03-N00
      5 1875  985 (N15W58 L=105)   466    119        1 (1846/ 965-1899/1013)   0.0 W54-W61 N13-N16
      6 1785  816 (N02W55 L=101)   940    168        1 (1757/ 739-1811/ 883) -50.2 W51-W59 N00-N06
      7  419  831 (S32E35 L= 11)  1094    318        3 ( 291/ 805- 538/ 863)  10.8 E46-E25 S35-S29
      8 1197  706 (S20W19 L= 65)   258     96        0 (1145/ 686-1248/ 718)  15.9 W15-W22 S22-S19
      9 1444  657 (S16W35 L= 81)   549    188        0 (1352/ 634-1525/ 684)  28.6 W28-W39 S18-S12
     10 1696  524 (S15W61 L=107)  2442    312        3 (1608/ 458-1762/ 611)  57.5 W53-W68 S19-S08
     11 1155  546 (S31W22 L= 68)   570    132        0 (1082/ 525-1225/ 573)  36.9 W17-W25 S33-S27
  

    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