Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1823 (N47E23 L= 86)   443     60        1 ( 949/1787- 992/1852) -29.7 E26-E19 N43-N49
      2 1122 1777 (N46E08 L=100)  1586    163        3 (1069/1691-1155/1847) -63.4 E13-E05 N38-N52
      3  776 1595 (N25E26 L= 82)   817    226        1 ( 689/1517- 848/1663) -26.6 E35-E19 N20-N28
      4  225 1146 (S09E58 L= 51)  4641    436       16 ( 131/1017- 310/1309)  57.1 E70-E48 S16-N01
      5 1353 1239 (N12W16 L=125)   295     86        0 (1318/1208-1396/1266)  51.3 W14-W18 N09-N14
      6  603 1077 (S09E26 L= 82)   311    136        0 ( 579/1018- 639/1128) -50.2 E28-E23 S12-S06
      7 1684 1067 (N08W41 L=150)   252     76        1 (1659/1041-1722/1089)   0.0 W39-W44 N06-N09
      8  396 1034 (S14E41 L= 68)   304     78        3 ( 373/1005- 418/1074)  76.0 E42-E39 S15-S11
      9  251  902 (S23E55 L= 54)  1647    276        3 ( 210/ 765- 290/1013)  68.2 E59-E50 S31-S16
     10 1015  976 (S09W00 L=109)   223     66        0 ( 974/ 961-1047/ 993)  33.7 E02-W02 S10-S07
     11 1815  888 (N01W57 L=165)   436    148        0 (1743/ 866-1868/ 906)   8.7 W49-W62 N00-N02
     12  378  761 (S31E43 L= 66)  1198    255        2 ( 360/ 644- 396/ 888)  76.0 E46-E40 S38-S22
     13 1750  572 (S18W64 L=173)  1327    122        2 (1706/ 527-1792/ 620)  66.8 W60-W69 S21-S14
     14 1742  491 (S23W70 L=179)   516     57        0 (1718/ 468-1768/ 516)   0.0 W65-W77 S23-S21
     15 1190  464 (S37W23 L=132)   302     92        1 (1176/ 420-1206/ 507) -50.2 W21-W26 S40-S34
     16 1559  368 (S35W61 L=169)  2131    274       11 (1489/ 300-1643/ 422)  55.0 W54-W69 S39-S29
  

    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