Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  994 1622 (N35E17 L=179)  3047    638        8 ( 752/1500-1233/1756)  -7.0 E39-W02 N30-N39
      2  399 1642 (N23E63 L=133)   581    201        1 ( 383/1591- 415/1676) -32.0 E66-E58 N20-N25
      3 1305 1493 (N33W08 L=204)   282     64        0 (1274/1476-1342/1510)  36.9 W06-W10 N31-N34
      4  627 1447 (N15E34 L=162)   401     69        0 ( 601/1410- 660/1490)  80.5 E35-E33 N12-N18
      5  391 1269 (N00E46 L=150)   552     87        1 ( 349/1241- 431/1302)  36.9 E49-E43 S01-N02
      6 1605 1072 (N13W36 L=232)   317     58        1 (1581/1035-1624/1111) -63.4 W33-W37 N11-N15
      7  760  994 (S09E15 L=181)   528    158        2 ( 732/ 928- 814/1065) -40.6 E17-E10 S12-S04
      8 1868 1001 (N15W61 L=257)   731    163        2 (1829/ 963-1898/1035) -12.5 W56-W66 N13-N16
      9  936  908 (S10E03 L=193)   343     93        1 ( 892/ 880- 990/ 933)  38.7 E05-W00 S11-S07
     10 1678  674 (S08W52 L=248)   578    118        0 (1663/ 627-1695/ 724) -39.8 W49-W55 S11-S05
     11 1341  620 (S19W29 L=225)   249     54        0 (1316/ 607-1377/ 636)  26.6 W26-W31 S19-S17
     12 1020  564 (S29W10 L=206)   408     88        3 ( 986/ 543-1064/ 596)   0.0 W07-W13 S30-S27
  

    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