Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1174 1856 (N65E08 L= 82)  3898    442        5 ( 983/1781-1351/1918)  -2.2 E39-W17 N58-N69
      2 1135 1680 (N45E04 L= 85)  1126    199        3 (1051/1623-1206/1722) -12.1 E12-W02 N41-N47
      3  551 1693 (N35E54 L= 35)   478     92        1 ( 507/1666- 595/1709)  33.7 E57-E50 N32-N37
      4  346 1507 (N19E62 L= 28)   458     68        2 ( 327/1482- 374/1541) -18.4 E65-E56 N17-N21
      5  308 1398 (N12E59 L= 30)   439     79        1 ( 292/1362- 328/1432) -45.0 E62-E57 N09-N13
      6 1514 1201 (N18W29 L=119)  1030    297        2 (1437/1079-1595/1320) -47.3 W22-W36 N11-N24
      7 1858 1074 (N15W63 L=152)  1006    139        2 (1833/1015-1886/1147) -63.4 W59-W67 N12-N20
      8  157  864 (S23E70 L= 20)   594     52        2 ( 141/ 840- 173/ 888) -45.0 E72-E66 S24-S21
      9 1211  735 (S15W16 L=106)   223     61        0 (1184/ 716-1244/ 755)  53.1 W14-W17 S16-S12
     10  557  706 (S27E27 L= 62)   835    173        2 ( 467/ 676- 635/ 729)  17.1 E34-E21 S28-S24
     11 1049  595 (S26W09 L= 98)   329    159        0 ( 990/ 550-1108/ 649)  56.3 W05-W11 S30-S21
     12  677  213 (S69E26 L= 63)   844     90        0 ( 651/ 177- 707/ 253)  79.7 E32-E21 S75-S63
  

    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