Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1629 (N37W18 L=239)   363     55        0 (1087/1601-1148/1655)  18.4 W15-W21 N34-N38
      2  361 1542 (N45E49 L=172)   397     41        0 ( 332/1530- 392/1555) -12.5 E53-E44 N44-N46
      3  671 1382 (N28E17 L=203)   508    106        1 ( 640/1336- 708/1421)  45.0 E19-E14 N25-N30
      4  424 1318 (N28E37 L=184)  2617    498        3 ( 210/1258- 617/1385)  -5.6 E61-E20 N23-N33
      5 1724  930 (S16W45 L=266)  1169    297        2 (1626/ 835-1823/1040)  25.3 W36-W55 S20-S11
      6 1273  813 (S16W12 L=232)   288     81        1 (1240/ 781-1313/ 835)  21.8 W09-W14 S17-S14
      7  735  768 (S10E22 L=199)   210     66        0 ( 713/ 735- 764/ 799)  36.9 E23-E19 S11-S08
      8 1462  732 (S24W23 L=244)   350     94        0 (1443/ 686-1481/ 786) -81.9 W22-W24 S27-S20
      9 1569  704 (S27W31 L=252)   495    107        0 (1519/ 677-1618/ 740)  14.0 W27-W35 S28-S25
     10 1251  682 (S23W08 L=229)   296    101        0 (1211/ 645-1300/ 703)   9.5 W05-W11 S25-S22
     11  414  544 (S18E52 L=169)  1515    389        2 ( 314/ 444- 501/ 668)  17.9 E70-E39 S23-S11
  

    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