Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  989 1580 (N42W13 L=255)   542     94        3 ( 948/1544-1022/1617)  18.4 W09-W15 N38-N43
      2 1699 1531 (N18W62 L=304)   613     98        1 (1663/1489-1728/1574) -63.4 W60-W65 N15-N21
      3 1494 1330 (N13W37 L=278)   233     53        0 (1467/1318-1527/1348) -45.0 W34-W38 N12-N15
      4  545 1175 (N25E28 L=214)  2323    524        5 ( 340/1089- 757/1265) -28.1 E44-E14 N16-N33
      5  272 1137 (N27E51 L=191)   388     70        1 ( 257/1108- 296/1169) -90.0 E52-E48 N25-N29
      6 1840  941 (S19W54 L=296)  2250    326        4 (1775/ 823-1881/1083)  40.9 W46-W62 S24-S11
      7  214  892 (N13E61 L=181)   724    193        3 ( 165/ 841- 259/ 941) -41.2 E65-E56 N09-N16
      8  993  796 (S05E07 L=235)   307    122        0 ( 930/ 767-1055/ 815)  -7.1 E10-E02 S06-S04
      9 1209  739 (S13W05 L=246)   196     58        0 (1196/ 709-1225/ 773)  76.0 W03-W05 S15-S11
     10 1398  647 (S24W15 L=256)   616    130        1 (1359/ 591-1426/ 703) -90.0 W12-W16 S27-S19
     11  292  564 (S08E65 L=177)   966    155        1 ( 246/ 488- 339/ 630) -68.2 E67-E62 S13-S03
     12 1671  550 (S40W37 L=278)   421     67        0 (1643/ 532-1699/ 564) -18.4 W33-W39 S41-S37
  

    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