Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1300 1619 (N36E00 L=112)  1693    288        0 (1187/1501-1405/1723) -24.2 E12-W09 N31-N41
      2  542 1664 (N19E51 L= 62)   567    130        1 ( 494/1640- 586/1713)  63.4 E54-E48 N16-N22
      3  584 1578 (N14E43 L= 69)   719    150        1 ( 518/1539- 664/1628)  45.0 E47-E39 N10-N18
      4  882 1489 (N16E21 L= 92)   219     68        0 ( 840/1479- 918/1501)  21.8 E23-E18 N14-N16
      5 1382 1427 (N26W10 L=123)   273     52        1 (1356/1402-1404/1453)  90.0 W08-W11 N24-N28
      6  240 1157 (S17E56 L= 57)   412     84        2 ( 224/1120- 262/1201) -78.7 E57-E54 S19-S14
      7 1292  618 (S22W28 L=140)   598    163        0 (1215/ 596-1372/ 642)  21.8 W22-W33 S23-S19
      8 1675  407 (S19W68 L=180)  2996    365        6 (1597/ 343-1745/ 492)  69.0 W62-W73 S25-S12
      9 1476  410 (S27W50 L=163)   415     31        1 (1450/ 394-1502/ 428)  45.0 W48-W52 S28-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