Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1229 1898 (N60W65 L= 31)   837    105        1 (1174/1878-1272/1921)  -2.0 W52-W82 N58-N61
      2  637 1701 (N59E24 L=302)   367     51        0 ( 614/1683- 662/1726)  -9.5 E27-E20 N57-N61
      3 1158 1674 (N45W27 L=352)   945    182        3 (1079/1634-1250/1700)  -3.6 W18-W34 N43-N46
      4  541 1673 (N58E37 L=289)   527     56        2 ( 515/1645- 566/1698)  63.4 E41-E33 N55-N60
      5 1301 1644 (N39W36 L=  2)   327     52        0 (1266/1623-1329/1663) -53.1 W33-W38 N37-N41
      6 1556 1617 (N31W56 L= 22)   712     64        2 (1532/1581-1591/1649)  45.0 W51-W60 N29-N33
      7 1373 1508 (N29W34 L=360)   290     93        0 (1329/1483-1410/1529) -51.3 W31-W36 N26-N31
      8  865 1388 (N31E04 L=322)  1988    536        6 ( 628/1310-1093/1480) -31.0 E21-W09 N22-N40
      9 1720 1444 (N17W60 L= 26)   443     96        1 (1671/1428-1761/1459) -20.6 W55-W64 N15-N18
     10 1784 1319 (N09W60 L= 26)   373     83        1 (1754/1299-1816/1347) -26.6 W57-W63 N06-N10
     11 1474 1269 (N13W33 L=358)   282     81        0 (1442/1232-1495/1308) -78.7 W31-W33 N10-N15
     12 1116 1089 (N09W06 L=332)   234     77        0 (1102/1044-1132/1127)  68.2 W04-W07 N06-N11
     13  441  995 (N13E38 L=288)  1092    300        1 ( 335/ 873- 534/1067) -42.9 E46-E32 N04-N18
     14 1427  871 (S09W22 L=348)   548    130        0 (1405/ 808-1448/ 938)  63.4 W20-W24 S13-S05
     15  919  619 (S16E13 L=312)   311    111        0 ( 878/ 587- 953/ 662)  33.7 E16-E10 S17-S13
     16  588  576 (S14E36 L=290)   263     72        0 ( 566/ 547- 619/ 611) -68.2 E36-E34 S16-S11
     17 1089  446 (S31E06 L=319)  1260    326        2 ( 996/ 378-1198/ 577)  30.5 E14-W03 S34-S23
     18  676  220 (S43E52 L=274)   916     99        0 ( 628/ 187- 718/ 249)  10.8 E65-E44 S45-S40
  

    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