Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  883 1831 (N64E20 L=  7)   640     85        0 ( 845/1810- 928/1849)  10.3 E25-E13 N61-N66
      2  660 1484 (N33E28 L=359)   280     66        2 ( 635/1464- 692/1511)  -9.5 E30-E24 N31-N34
      3 1337 1406 (N28W22 L= 49)   341    100        0 (1294/1371-1379/1437) -38.7 W19-W25 N25-N29
      4  980 1332 (N23E03 L= 23)   329    150        0 ( 917/1280-1028/1370) -40.6 E07-W00 N19-N25
      5 1810  845 (S09W59 L= 85)   614     96        0 (1797/ 794-1822/ 903)  90.0 W57-W60 S12-S05
      6 1753  856 (S08W52 L= 79)   570     45        1 (1728/ 829-1775/ 886) -63.4 W49-W54 S10-S06
      7 1544  750 (S14W35 L= 61)  1153    362        1 (1391/ 703-1645/ 812)  15.5 W24-W43 S16-S10
      8 1025  694 (S17E00 L= 26)   851    190        1 ( 934/ 652-1112/ 729)  21.8 E05-W05 S19-S14
      9 1525  606 (S24W36 L= 62)  1298    302        0 (1436/ 541-1571/ 709)  59.0 W29-W39 S28-S16
     10  928  602 (S23E07 L= 20)   246     79        0 ( 916/ 557- 939/ 637) -90.0 E07-E05 S26-S20
     11  504  496 (S32E41 L=345)  2638    373        8 ( 420/ 420- 605/ 599) -45.0 E47-E33 S37-S24
  

    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