Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1819 (N42E30 L=174)  1053    212        0 ( 937/1758-1025/1881) -31.0 E38-E23 N37-N47
      2  420 1645 (N16E60 L=145)   551    152        0 ( 405/1598- 433/1685) -45.0 E64-E56 N12-N18
      3  919 1412 (N12E16 L=189)  1188    309        0 ( 874/1260- 982/1535) -81.0 E19-E12 N02-N21
      4 1805 1199 (N26W47 L=252)   350     34        0 (1790/1176-1822/1228) -33.7 W45-W48 N24-N27
      5 1100 1177 (N03W00 L=205)   980    282        1 ( 950/1098-1214/1214)   3.6 E08-W08 N00-N05
      6 1025  903 (S14W03 L=208)   878    216        4 ( 951/ 861-1144/ 949)  45.0 W00-W08 S17-S08
      7 1745  892 (N06W48 L=253)   449    103        0 (1708/ 870-1790/ 914)   7.1 W45-W53 N04-N07
      8  769  591 (S39E04 L=201)  1991    490        4 ( 574/ 548-1020/ 628)  17.4 E19-W13 S43-S32
      9 1355  548 (S24W35 L=239)  1346    368        0 (1190/ 519-1520/ 586)  24.4 W24-W46 S29-S19
     10 1484  473 (S24W47 L=252)   305     56        1 (1460/ 448-1512/ 495)  71.6 W45-W47 S25-S22
     11 1196  281 (S45W41 L=246)   409     70        1 (1163/ 259-1224/ 302)  11.3 W38-W43 S46-S43
  

    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