Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1779 (N32E49 L=230)   569    136        0 ( 683/1723- 768/1828) -68.2 E52-E46 N26-N36
      2 1070 1701 (N36E20 L=258)  2052    412        2 ( 932/1558-1153/1770) -22.2 E31-E09 N28-N41
      3  536 1697 (N21E57 L=222)  1131    164        3 ( 478/1657- 592/1762)  83.7 E62-E54 N17-N26
      4  595 1672 (N21E50 L=229)   381     55        1 ( 570/1653- 625/1692)   0.0 E52-E46 N19-N22
      5  734 1618 (N21E37 L=242)   439    114        0 ( 718/1548- 751/1664) -49.4 E39-E33 N16-N23
      6  319 1268 (S09E53 L=226)   829    168        5 ( 248/1232- 351/1329)  40.6 E59-E49 S12-S05
      7  172 1143 (S19E67 L=212)  4474    409       10 (  94/ 959- 220/1296)  57.5 E84-E61 S31-S09
      8  274 1176 (S16E55 L=224)   608     57        2 ( 249/1142- 303/1212) -45.0 E57-E51 S17-S13
      9 1636  983 (N10W37 L=316)   847    331        0 (1598/ 852-1662/1144) -53.7 W30-W41 N03-N18
     10  313  765 (S40E46 L=233)   866    127        2 ( 263/ 722- 375/ 792)  21.8 E53-E38 S43-S37
     11 1333  735 (S12W26 L=305)   334    109        0 (1290/ 700-1380/ 761)  59.0 W24-W28 S15-S09
     12  799  589 (S37E00 L=278)   226      6        0 ( 780/ 572- 819/ 606)   0.0 E01-W01 S38-S35
  

    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