Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1286 1682 (N35W36 L=177)  1720    348        5 (1185/1615-1409/1716) -14.0 W27-W44 N31-N37
      2  710 1646 (N49E07 L=134)   273     75        0 ( 684/1626- 740/1672) -38.7 E09-E04 N47-N51
      3  572 1546 (N45E21 L=120)   392    122        2 ( 523/1525- 634/1563)  -5.2 E25-E14 N43-N46
      4 1488 1038 (S07W27 L=168)  1275    348        3 (1415/ 870-1556/1200)  41.2 W19-W35 S14-N00
      5 1460  752 (S23W21 L=162)   726    209        2 (1344/ 707-1530/ 817)   8.1 W12-W26 S24-S20
      6 1082  506 (S29E09 L=132)  1283    338        3 (1014/ 397-1111/ 668)  50.9 E16-E03 S35-S19
      7  664  448 (S24E37 L=104)  1214    412        2 ( 636/ 301- 707/ 575)  46.7 E45-E29 S33-S16
      8 1493  522 (S39W21 L=161)   680    123        1 (1441/ 485-1542/ 559)  14.0 W16-W24 S40-S36
  

    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