Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  450 1577 (N21E54 L= 42)   458    113        1 ( 414/1557- 492/1590)  18.4 E57-E50 N18-N21
      2  421 1276 (N02E44 L= 51)   284     66        0 ( 390/1262- 456/1291)  21.8 E46-E41 N00-N03
      3 1216 1208 (N13W08 L=103)   501    168        1 (1154/1157-1284/1291) -26.6 W02-W12 N11-N16
      4  598  983 (S12E25 L= 70)   327    132        0 ( 581/ 922- 614/1055) -66.8 E27-E24 S15-S08
      5 1786  928 (N08W54 L=149)   403    104        0 (1752/ 906-1816/ 949)  -8.1 W50-W57 N07-N09
      6 1374  885 (S03W24 L=119)   220     60        0 (1353/ 853-1394/ 920)  76.0 W22-W24 S04-N00
      7  897  660 (S26E00 L= 95)  1188    543        0 ( 657/ 557-1147/ 750)  39.4 E15-W13 S38-S15
      8 1393  624 (S18W32 L=127)   235     80        0 (1353/ 594-1429/ 644)  45.0 W30-W34 S20-S16
      9 1546  564 (S18W45 L=141)   407     77        1 (1520/ 535-1570/ 599)  76.0 W43-W47 S19-S15
     10 1576  460 (S24W53 L=148)   371    117        1 (1533/ 445-1621/ 471)  23.2 W49-W56 S24-S21
     11 1408  269 (S40W53 L=148)  1713    292        1 (1288/ 247-1523/ 301)  33.7 W44-W61 S45-S35
  

    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