Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  592 1806 (N31E57 L=123)   698    124        3 ( 563/1787- 622/1833)   0.0 E60-E53 N28-N33
      2 1262 1686 (N40E05 L=175)   533     73        0 (1228/1651-1292/1714) -23.2 E08-E01 N38-N42
      3  632 1651 (N20E43 L=137)   318     50        0 ( 603/1638- 662/1665)  11.3 E45-E40 N19-N21
      4  699 1534 (N14E33 L=147)   465    159        2 ( 676/1471- 730/1600) -53.1 E36-E29 N10-N18
      5  390 1019 (S22E40 L=140)   494    183        1 ( 323/ 980- 461/1051)  37.9 E45-E34 S25-S18
      6 1454 1011 (N04W25 L=204)   508    167        3 (1379/ 988-1535/1031)  24.0 W20-W29 N01-N06
      7  699  839 (S26E15 L=165)   414    181        0 ( 612/ 808- 781/ 874)  35.0 E20-E10 S29-S21
      8 1562  745 (S07W40 L=220)   924    226        3 (1515/ 636-1606/ 824) -80.5 W36-W42 S14-S02
      9 1325  698 (S17W27 L=206)  1417    363        4 (1158/ 645-1481/ 773)  41.2 W17-W34 S22-S08
     10 1010  460 (S40W18 L=197)  1621    380        6 ( 853/ 430-1182/ 495)  19.2 W06-W29 S44-S35
     11 1511  400 (S27W53 L=233)  1629    390        1 (1352/ 370-1651/ 430)  29.7 W41-W63 S32-S20
  

    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