Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1277 1176 (N06W18 L=146)   197     73        0 (1249/1139-1307/1214)  36.9 W15-W20 N04-N07
      2 1414 1147 (N01W25 L=153)   293     84        1 (1375/1103-1451/1180) -53.1 W23-W26 S02-N03
      3  523 1115 (N21E28 L=100)   283     83        0 ( 500/1073- 557/1156)  31.0 E30-E24 N19-N23
      4 1554 1013 (S10W31 L=159)   321    146        1 (1503/1000-1623/1027) -32.0 W27-W35 S12-S07
      5 1203  967 (S03W08 L=137)   185     82        2 (1171/ 922-1234/ 998) -80.5 W07-W09 S06-N00
      6  407  910 (N12E41 L= 88)   317    162        0 ( 345/ 892- 466/ 934) -11.3 E45-E35 N11-N13
      7  165  876 (N15E64 L= 64)   745    157        3 ( 131/ 813- 193/ 918) -60.3 E67-E62 N10-N17
      8 1775  846 (S26W45 L=173)   415    129        1 (1747/ 790-1802/ 913)  35.5 W41-W48 S29-S23
      9  856  647 (S13E19 L=109)   236     67        0 ( 825/ 623- 885/ 671)  21.8 E21-E16 S13-S11
     10 1349  570 (S31W08 L=136)   509     60        1 (1304/ 537-1387/ 606)   9.5 W04-W10 S32-S28
     11 1066  433 (S31E14 L=114)   297     66        0 (1047/ 402-1083/ 463)  53.1 E16-E12 S33-S29
     12  822  378 (S28E32 L= 96)   351    139        0 ( 807/ 321- 833/ 430)  49.4 E35-E28 S32-S25
  

    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