Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  807 1626 (N23E32 L=206)   375    121        2 ( 778/1579- 845/1671)  90.0 E33-E30 N19-N27
      2  526 1613 (N16E50 L=188)   350    100        1 ( 499/1584- 552/1642) -15.9 E53-E46 N14-N17
      3  452 1483 (N06E50 L=189)   393    109        0 ( 412/1467- 493/1497)   0.0 E53-E46 N05-N07
      4 1366 1347 (N22W11 L=249)   432     65        1 (1338/1313-1395/1383) -36.9 W09-W13 N20-N24
      5  300 1172 (S15E51 L=187)  1196    311        0 ( 201/1085- 365/1298)  61.9 E58-E46 S21-S06
      6  206 1203 (S15E62 L=177)   814    161        1 ( 167/1152- 253/1258)  58.0 E65-E58 S18-S10
      7  360  902 (S30E40 L=198)   345     75        0 ( 343/ 854- 378/ 942) -59.0 E42-E38 S32-S27
      8 1055  870 (S13W06 L=244)   319     76        2 (1026/ 834-1093/ 904)  68.2 W04-W07 S15-S10
      9  697  819 (S26E14 L=225)   306     92        3 ( 674/ 766- 719/ 862) -59.0 E15-E11 S28-S23
     10  194  791 (S40E56 L=182)   681    146        0 ( 164/ 771- 233/ 826)  -7.1 E59-E51 S41-S38
     11 1009  668 (S26W10 L=248)   238     61        0 ( 984/ 638-1037/ 706) -21.8 W07-W12 S26-S24
     12 1430  666 (S13W34 L=273)   414     99        0 (1387/ 623-1483/ 705)  66.8 W32-W36 S16-S09
     13  832  652 (S32E01 L=238)   272     50        1 ( 797/ 626- 861/ 670) -14.0 E03-W01 S32-S30
     14 1708  482 (S13W67 L=305)   653    126        1 (1683/ 447-1724/ 523) -31.0 W63-W69 S14-S11
     15  847  389 (S49W11 L=250)   286     63        0 ( 830/ 362- 860/ 426) -71.6 W09-W13 S51-S45
     16  975  288 (S51W28 L=266)  1553    326        5 ( 923/ 220-1044/ 380) -30.5 W20-W37 S55-S44
     17 1273  335 (S37W42 L=281)   761    109        3 (1235/ 306-1317/ 367)  45.0 W38-W45 S39-S33
  

    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