Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  547 1702 (N46E48 L=129)  3693    552       10 ( 422/1625- 686/1761)  21.8 E57-E34 N39-N50
      2  854 1629 (N42E17 L=160)   355     73        0 ( 822/1603- 882/1660)  45.0 E19-E15 N40-N44
      3 1257 1545 (N37W16 L=193)   250     41        1 (1235/1527-1283/1565) -33.7 W14-W17 N36-N38
      4  257 1377 (N21E61 L=115)   695     99        2 ( 236/1329- 280/1430) -81.9 E63-E59 N17-N24
      5  444 1360 (N21E43 L=134)   573     95        2 ( 419/1325- 474/1406) -51.3 E45-E39 N18-N24
      6  390 1249 (N14E45 L=132)  1025    170        3 ( 340/1187- 434/1295) -56.3 E49-E41 N10-N16
      7  203 1188 (N09E64 L=113)   486    128        1 ( 186/1140- 224/1236) -50.2 E66-E59 N05-N11
      8  375 1162 (N08E45 L=132)   883    306        1 ( 303/1130- 464/1212) -11.3 E52-E37 N06-N10
      9 1875  849 (S07W68 L=244)   753    115        3 (1853/ 805-1896/ 885) -59.0 W63-W71 S09-S04
     10  385  805 (S14E43 L=133)  1290    234        3 ( 313/ 753- 452/ 875) -28.6 E49-E38 S16-S09
     11 1442  582 (S24W31 L=208)   267     88        0 (1403/ 567-1478/ 602) -18.4 W27-W33 S25-S22
     12 1251  431 (S36W20 L=197)   927    187        1 (1179/ 406-1346/ 468)  -4.4 W13-W26 S37-S33
     13  861  402 (S40E10 L=166)  1343    341        3 ( 708/ 353-1008/ 441)  14.0 E23-W01 S44-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