Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  931 1569 (N31E18 L= 67)   552    256        1 ( 840/1509-1029/1642) -16.4 E27-E10 N27-N34
      2 1246 1611 (N41W03 L= 88)   391    102        0 (1206/1595-1298/1626)  18.4 W00-W06 N38-N42
      3  491 1555 (N21E50 L= 35)  1647    277        8 ( 373/1509- 591/1604)  22.4 E59-E42 N16-N25
      4 1596 1377 (N32W34 L=119)   422     62        1 (1570/1350-1622/1411) -45.0 W31-W36 N30-N34
      5 1333 1268 (N19W14 L=100)   339     96        3 (1295/1243-1377/1292)  -9.5 W11-W17 N17-N20
      6 1389 1161 (N14W20 L=105)   238     90        0 (1357/1119-1422/1181) -31.0 W17-W22 N11-N15
      7  192  990 (S18E62 L= 23)  2902    336        6 ( 114/ 862- 255/1128)  33.7 E83-E55 S27-S09
      8 1781 1043 (N15W52 L=138)   391     73        0 (1762/1011-1798/1080)  68.2 W50-W54 N12-N17
      9 1049 1014 (S01W02 L= 87)   644    186        1 ( 958/ 984-1120/1042)  26.6 E03-W06 S04-N00
     10 1123  844 (S10W09 L= 95)   195    119        0 (1071/ 830-1171/ 863)   8.1 W05-W12 S11-S08
     11  783  829 (S17E11 L= 74)   661    293        3 ( 657/ 792- 917/ 857)  28.1 E18-E03 S22-S13
     12 1398  690 (S14W30 L=116)   923    296        2 (1278/ 602-1498/ 765)  54.5 W25-W35 S21-S07
     13  393  652 (S37E41 L= 45)  1372    335        2 ( 305/ 552- 508/ 743)  35.8 E55-E29 S46-S28
     14 1712  643 (S10W57 L=143)  1781    257        3 (1604/ 554-1785/ 719)  45.0 W49-W64 S16-S04
     15 1737  604 (S12W62 L=148)   521     90        1 (1711/ 572-1772/ 643)  63.4 W59-W64 S14-S08
     16 1468  506 (S24W43 L=128)   802    286        2 (1398/ 467-1526/ 579)  60.9 W37-W46 S26-S17
     17 1180  526 (S28W21 L=107)   225     33        0 (1156/ 511-1208/ 540)  14.0 W19-W23 S29-S27
     18  959  475 (S37W08 L= 94)   271     78        1 ( 922/ 459-1002/ 497)  38.7 W05-W10 S38-S34
  

    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