Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1714 1522 (N12W64 L=131)  1758    176        4 (1661/1455-1770/1593) -84.3 W60-W67 N07-N17
      2  291 1475 (N52E55 L= 11)  1958    200        3 ( 197/1441- 380/1500)  -2.7 E82-E40 N49-N54
      3 1051 1461 (N29W14 L= 81)  1110    297        2 ( 898/1419-1185/1497) -30.5 W05-W22 N24-N34
      4 1450 1427 (N15W37 L=104)   578    260        0 (1334/1401-1570/1444) -20.6 W30-W46 N12-N19
      5  836 1377 (N30E02 L= 65)   217     48        1 ( 809/1359- 864/1397) -56.3 E03-W00 N28-N31
      6  457 1170 (N28E33 L= 34)   264     69        0 ( 415/1158- 495/1179) -15.9 E36-E29 N26-N28
      7  388  984 (N18E41 L= 25)  1456    391        2 ( 357/ 840- 448/1171)  81.9 E44-E37 N08-N29
      8 1812  927 (S23W51 L=117)   526     73        1 (1785/ 889-1836/ 963)  45.0 W47-W53 S25-S21
      9 1219  654 (S21W01 L= 67)   462    227        0 (1102/ 642-1322/ 668) -26.6 E05-W07 S23-S17
     10  425  447 (S14E59 L=  7)  1021    185        0 ( 401/ 363- 467/ 515)  25.3 E69-E50 S19-S10
     11  646  351 (S25E47 L= 19)   563     93        1 ( 616/ 315- 677/ 402)  29.1 E51-E42 S27-S21
     12  907  334 (S33E30 L= 36)   436     82        1 ( 881/ 305- 930/ 376)  29.7 E33-E26 S35-S30
  

    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