Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  849 1539 (N36E15 L=225)   267     84        0 ( 814/1522- 900/1552)  26.6 E17-E11 N34-N37
      2  336 1504 (N31E62 L=179)   579     72        1 ( 316/1474- 366/1536) -15.3 E67-E56 N28-N32
      3 1601 1355 (N25W42 L=282)   269     82        0 (1555/1339-1636/1374) -15.9 W37-W44 N23-N25
      4  412 1265 (N16E44 L=197)   635    185        4 ( 376/1191- 464/1311)  74.1 E46-E40 N11-N18
      5  297  723 (S19E54 L=186)   372     84        0 ( 265/ 709- 335/ 736)   0.0 E58-E50 S20-S18
      6  324  645 (S24E54 L=187)   766    167        2 ( 302/ 570- 369/ 705)  45.0 E59-E47 S29-S19
      7 1468  396 (S39W40 L=281)   510    124        0 (1408/ 372-1518/ 423)  26.6 W34-W44 S40-S36
      8  981  362 (S42E01 L=239)   782    136        2 ( 937/ 322-1019/ 396)  45.0 E05-W01 S46-S39
      9  696  359 (S44E27 L=214)   307     64        0 ( 667/ 344- 725/ 379) -21.8 E29-E24 S45-S42
     10  796  330 (S46E18 L=222)   390     43        2 ( 771/ 302- 819/ 352) -90.0 E20-E15 S48-S44
  

    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