Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  863 1512 (N19E24 L= 26)   231     68        0 ( 828/1498- 898/1526)  11.3 E25-E20 N17-N19
      2 1424 1426 (N30W13 L= 63)   375    152        0 (1379/1373-1470/1452)  60.3 W10-W16 N25-N33
      3 1856 1248 (N34W58 L=108)   465     43        1 (1837/1219-1880/1267)  45.0 W54-W62 N31-N36
      4 1819 1025 (N19W54 L=104)  3362    433        9 (1749/ 898-1880/1193) -55.5 W47-W62 N11-N28
      5 1618 1137 (N18W33 L= 83)   259     52        0 (1590/1115-1644/1155)   0.0 W30-W35 N17-N19
      6 1579 1045 (N12W32 L= 82)   346     70        1 (1549/1007-1599/1083) -53.1 W30-W34 N10-N14
      7  129 1025 (S26E73 L=337)   912     83        1 ( 109/ 978- 147/1056) -45.0 E78-E69 S28-S23
      8  252  915 (S30E53 L=357)   636    163        0 ( 199/ 876- 302/ 945)  26.6 E60-E47 S33-S27
      9  489  548 (S49E27 L= 23)   348     25        1 ( 464/ 529- 510/ 572)   0.0 E29-E24 S50-S46
     10  826  503 (S41W04 L= 54)   252     67        0 ( 797/ 484- 862/ 528)   0.0 W01-W06 S41-S39
     11 1038  495 (S34W18 L= 68)   510     72        0 ( 998/ 479-1088/ 518)  33.7 W15-W21 S36-S32
     12 1370  338 (S33W50 L=101)   525     84        1 (1324/ 312-1399/ 365)   5.7 W44-W54 S34-S31
     13 1148  171 (S53W57 L=107)   531     60        1 (1123/ 142-1167/ 205) -21.0 W51-W64 S54-S49
  

    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