Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  468 1444 (N46E33 L= 56)   526    134        1 ( 405/1423- 534/1471)  -8.1 E40-E26 N44-N46
      2  981 1031 (N09E02 L= 87)   251     85        1 ( 960/ 999- 997/1077)  45.0 E04-W00 N06-N10
      3  272  913 (N15E55 L= 35)  1326    230        2 ( 224/ 832- 319/1024) -61.7 E58-E50 N09-N22
      4  203  867 (N13E64 L= 26)   751    119        2 ( 182/ 814- 226/ 932) -90.0 E67-E61 N09-N17
      5 1379  701 (S20W14 L=103)  1691    519        2 (1210/ 537-1531/ 880)  25.7 W00-W27 S27-S13
      6 1627  793 (S22W33 L=122)   450    180        0 (1607/ 711-1648/ 880)  63.4 W30-W35 S26-S16
      7 1093  585 (S20E06 L= 83)  1394    517        4 (1035/ 438-1177/ 791)  56.3 E12-W03 S30-S09
      8 1706  645 (S35W40 L=129)   500     51        1 (1685/ 621-1729/ 668)  26.6 W37-W42 S36-S32
      9  266  612 (S03E66 L= 23)   640     72        0 ( 248/ 582- 294/ 648) -76.0 E68-E63 S05-N00
     10 1434  533 (S33W15 L=104)   747    209        1 (1373/ 465-1470/ 628)  45.0 W08-W17 S36-S27
     11  730  262 (S35E46 L= 43)  1983    251        2 ( 674/ 163- 776/ 381)  41.0 E63-E40 S44-S24
  

    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