Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  931 1558 (N39W11 L=163)   747    197        2 ( 872/1476-1016/1609)  21.0 W05-W18 N34-N40
      2 1304 1512 (N24W32 L=184)  1289    253        2 (1185/1468-1419/1560) -47.5 W25-W37 N18-N30
      3  286 1208 (N34E46 L=106)   310     86        1 ( 258/1186- 318/1227) -21.8 E48-E43 N32-N35
      4  450 1153 (N27E33 L=119)   550    152        1 ( 385/1132- 515/1176) -16.7 E37-E27 N25-N29
      5  686 1082 (N17E18 L=134)   631    138        1 ( 629/1062- 760/1114) -26.6 E21-E13 N15-N20
      6 1321 1052 (S01W17 L=169)   912    216        3 (1218/1011-1431/1086) -17.1 W10-W23 S04-N01
      7 1292  975 (S04W13 L=165)   304    109        0 (1243/ 946-1347/ 993)  -8.1 W09-W16 S05-S03
      8 1508  891 (S15W24 L=176)   509    112        1 (1460/ 864-1557/ 928)   0.0 W20-W28 S16-S13
      9 1222  732 (S16W03 L=155)   738    201        0 (1190/ 655-1267/ 834)  74.7 W01-W06 S21-S10
     10  904  589 (S16E19 L=133)   654    151        1 ( 808/ 553- 965/ 620) -48.8 E23-E15 S19-S11
     11  584  432 (S17E45 L=107)   778    180        3 ( 535/ 399- 639/ 463)   0.0 E49-E40 S19-S15
     12  591  354 (S23E49 L=103)  1319    318        4 ( 481/ 278- 695/ 409) -60.3 E56-E45 S30-S16
     13  898  166 (S46E44 L=108)  1094    172        2 ( 874/ 129- 925/ 205)  35.0 E49-E39 S49-S41
  

    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