Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  551 1544 (N21E44 L= 42)   538    147        0 ( 500/1524- 607/1569)   0.0 E48-E39 N19-N22
      2  314 1529 (N15E65 L= 20)   502     81        0 ( 295/1495- 336/1565)  90.0 E67-E63 N13-N18
      3  406 1246 (S00E45 L= 41)   277     80        1 ( 385/1218- 437/1276) -21.8 E46-E41 S01-N01
      4 1345 1182 (N14W16 L=102)   442    158        1 (1295/1130-1412/1246) -29.1 W12-W21 N11-N17
      5 1860 1098 (N20W59 L=145)   418     92        0 (1831/1084-1884/1116)   8.1 W55-W62 N19-N21
      6  731  945 (S12E16 L= 70)   454    176        1 ( 680/ 872- 754/1014) -53.1 E20-E14 S16-S07
      7 1509  850 (S02W33 L=119)   480    108        1 (1479/ 802-1533/ 915)  90.0 W31-W34 S05-N02
      8 1877  879 (N07W67 L=153)   838    143        1 (1844/ 843-1904/ 908) -18.4 W61-W72 N05-N08
      9 1371  822 (S06W25 L=111)   219     58        1 (1356/ 793-1383/ 853) -56.3 W23-W26 S08-S04
     10  972  603 (S28W06 L= 91)  1002    387        1 ( 807/ 554-1173/ 654)  33.7 E04-W17 S34-S20
     11 1612  525 (S19W53 L=138)   867    167        1 (1534/ 482-1669/ 575)  52.1 W47-W55 S23-S14
     12 1622  441 (S24W58 L=144)   420    106        2 (1580/ 422-1652/ 459)  26.6 W54-W61 S25-S22
     13 1462  254 (S40W58 L=144)   822    198        1 (1430/ 238-1504/ 274)   0.0 W54-W63 S41-S38
  

    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