Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1543 (N35W09 L=131)   836    243        1 ( 865/1510-1069/1567) -11.3 W01-W16 N32-N36
      2 1596 1490 (N16W50 L=172)   413     52        0 (1570/1469-1621/1512) -56.3 W46-W51 N14-N17
      3 1467 1454 (N17W38 L=160)   252     80        0 (1453/1412-1488/1492)  45.0 W36-W41 N14-N19
      4 1285 1322 (N14W22 L=144)   403    132        1 (1227/1306-1355/1342) -14.0 W18-W26 N12-N15
      5  288 1260 (N32E49 L= 73)  1879    273        3 ( 199/1158- 387/1335) -36.9 E61-E38 N23-N38
      6  738 1187 (N18E13 L=109)  1559    432        1 ( 553/1123- 918/1293)   2.2 E26-W00 N15-N20
      7  140 1163 (N29E69 L= 53)   724     58        2 ( 121/1133- 157/1193) -76.0 E72-E65 N26-N30
      8 1667  811 (S24W37 L=159)  1298    319        0 (1611/ 675-1723/ 946)  51.8 W31-W43 S31-S17
      9 1856  784 (S32W59 L=181)   807     86        1 (1822/ 754-1889/ 822)  18.4 W54-W63 S33-S29
     10  511  553 (S14E44 L= 78)   668    120        1 ( 491/ 479- 531/ 615)  60.9 E48-E41 S18-S09
     11 1483  583 (S34W20 L=142)   294     39        1 (1458/ 561-1507/ 602)   0.0 W18-W22 S36-S32
     12 1023  538 (S26E12 L=111)   261     40        1 (1000/ 513-1044/ 562)  45.0 E12-E09 S27-S24
  

    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