Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1515 1549 (N20W48 L=245)   771    198        2 (1443/1517-1592/1570) -18.4 W42-W54 N18-N23
      2  781 1389 (N34E05 L=193)   222     65        0 ( 755/1355- 800/1429)  36.9 E06-E02 N31-N35
      3  543 1307 (N35E23 L=174)   502    142        1 ( 486/1264- 615/1347)   5.2 E28-E17 N33-N35
      4  964 1185 (N17W01 L=198)   206     64        2 ( 930/1155- 997/1204) -63.4 W00-W02 N14-N18
      5  418 1118 (N27E36 L=161)   354     66        1 ( 385/1098- 447/1138)   0.0 E38-E33 N25-N27
      6  368  949 (N18E43 L=154)   293     67        1 ( 345/ 920- 395/ 982) -68.2 E44-E41 N15-N20
      7 1024  526 (S22E14 L=183)   242     72        1 ( 996/ 500-1060/ 556) -63.4 E15-E12 S24-S19
      8  493  442 (S14E52 L=146)  1633    307        5 ( 417/ 356- 542/ 529) -79.7 E56-E45 S20-S09
  

    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