Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1433 1689 (N33W48 L=187)   357     91        0 (1390/1661-1472/1712) -59.0 W45-W49 N30-N35
      2 1619 1640 (N26W61 L=200)   437     32        0 (1594/1620-1642/1660)   0.0 W57-W64 N24-N27
      3  196 1271 (N31E61 L= 78)   493     75        0 ( 167/1244- 228/1292)   7.1 E65-E56 N30-N32
      4  709  975 (N05E19 L=120)   200     55        0 ( 683/ 951- 738/ 995)  26.6 E21-E17 N03-N05
      5 1641  937 (S15W37 L=176)   441    104        0 (1599/ 896-1671/ 981)  29.7 W33-W40 S16-S12
      6 1802  898 (S21W52 L=191)   983    223        1 (1707/ 837-1876/ 969)  12.5 W42-W60 S22-S17
      7  511  637 (S12E40 L= 99)  1146    197        3 ( 482/ 545- 545/ 730)  48.0 E45-E35 S16-S06
      8 1716  654 (S35W44 L=183)   309     38        1 (1702/ 631-1735/ 679) -53.1 W42-W46 S36-S32
      9 1451  582 (S33W20 L=159)   629    150        1 (1389/ 542-1519/ 621) -38.7 W15-W25 S37-S29
     10 1615  557 (S39W35 L=173)   386     89        1 (1588/ 530-1651/ 596)  31.0 W32-W37 S40-S37
     11 1245  479 (S36W03 L=142)   778    238        1 (1141/ 409-1329/ 563)  21.8 E05-W10 S39-S31
     12 1104  485 (S32E06 L=133)   651    133        1 (1080/ 429-1132/ 539)  66.8 E08-E03 S36-S29
  

    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