Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1151 1772 (N50E13 L= 51)  4172    230        8 (1071/1614-1236/1916) -37.1 E34-E01 N38-N64
      2  384 1540 (N15E59 L=  6)  2062    231        6 ( 339/1441- 433/1647) -45.0 E69-E53 N08-N21
      3  609 1220 (N00E29 L= 36)  1265    253        2 ( 542/1120- 689/1323)  72.9 E31-E26 S06-N07
      4 1362 1269 (N19W15 L= 80)   243     99        0 (1310/1247-1400/1305)  -9.5 W11-W17 N18-N20
      5  926 1107 (N00E07 L= 57)   193     90        0 ( 915/1058- 936/1150) -63.4 E08-E06 S02-N02
      6  818 1000 (S08E11 L= 53)   356    130        0 ( 761/ 953- 868/1032)  56.3 E13-E08 S11-S05
      7 1781  781 (N01W57 L=122)  3374    300        4 (1715/ 660-1845/ 928) -25.6 W47-W70 S04-N07
      8 1435  859 (S03W28 L= 92)   549    116        0 (1392/ 820-1479/ 897)  59.0 W25-W30 S05-N00
      9  433  693 (S35E35 L= 30)   376     89        0 ( 389/ 657- 459/ 715) -18.4 E38-E32 S36-S33
     10 1176  172 (S54W51 L=116)   463     66        0 (1153/ 144-1197/ 201) -31.0 W46-W57 S56-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