Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  781 1712 (N47E28 L=253)   352     18        0 ( 757/1696- 808/1727)   0.0 E31-E25 N45-N48
      2 1116 1622 (N42W03 L=284)  2823    555        4 ( 869/1511-1337/1686) -11.9 E18-W20 N35-N45
      3 1656 1513 (N37W51 L=332)  2146    279        1 (1566/1439-1761/1554)   4.6 W39-W65 N31-N39
      4 1411 1513 (N36W26 L=307)   515     67        0 (1383/1477-1445/1538)  31.0 W23-W29 N33-N37
      5  950 1188 (N12E06 L=275)  1752    401        1 ( 915/ 991- 971/1373) -75.4 E09-E03 N00-N23
      6 1651 1293 (N22W43 L=324)   405     82        0 (1630/1251-1670/1329) -68.2 W41-W44 N19-N24
      7 1723 1286 (N22W50 L=331)   397     78        0 (1713/1239-1735/1327)  63.4 W48-W52 N18-N24
      8  300 1162 (N06E52 L=229)  1859    347        1 ( 266/1010- 326/1326) -74.5 E58-E48 S03-N15
      9  469 1226 (N11E38 L=243)   988    319        2 ( 423/1139- 572/1318)  63.4 E41-E31 N05-N17
     10  503  756 (S18E33 L=247)  1445    448        0 ( 455/ 561- 542/ 950)  74.9 E38-E30 S31-S05
     11  210  894 (S11E60 L=221)  1043    117        2 ( 156/ 840- 242/ 939) -29.1 E67-E56 S14-S08
     12  945  621 (S24E03 L=278)   476    123        0 ( 882/ 593-1001/ 655)  33.7 E06-W00 S25-S20
     13 1734  570 (S23W60 L=341)   970    107        2 (1688/ 541-1790/ 594)  -7.1 W53-W69 S25-S21
     14 1557  465 (S31W46 L=327)  1812    206        3 (1462/ 430-1652/ 503)  13.2 W37-W55 S33-S28
     15 1334  380 (S39W30 L=311)   307     42        1 (1303/ 366-1362/ 398)   0.0 W27-W32 S39-S37
     16 1113  314 (S46W13 L=294)   348     66        1 (1069/ 300-1142/ 333)   8.1 W09-W16 S47-S43
     17  631  297 (S52E34 L=247)   461     59        0 ( 606/ 272- 653/ 317) -90.0 E38-E30 S54-S49
  

    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