Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1446 1557 (N39W34 L=123)   487    123        1 (1405/1514-1499/1593) -45.0 W31-W37 N35-N41
      2  308 1532 (N33E70 L= 20)   709     89        1 ( 293/1486- 322/1576) -31.0 E76-E65 N29-N35
      3 1693 1456 (N31W55 L=144)   587    100        1 (1652/1419-1728/1489) -39.8 W51-W58 N28-N33
      4  581 1449 (N29E34 L= 55)   349    105        1 ( 554/1419- 625/1477)  36.9 E36-E31 N26-N30
      5 1249 1370 (N25W15 L=104)   272    102        0 (1211/1337-1289/1397)  33.7 W11-W17 N23-N27
      6 1483 1278 (N20W31 L=120)   532    184        0 (1410/1216-1543/1332) -41.6 W25-W34 N15-N23
      7 1862 1236 (N16W68 L=157)  1724    199        4 (1830/1154-1887/1312)  56.3 W62-W74 N11-N20
      8 1754  855 (S07W53 L=142)   609     26        1 (1727/ 831-1780/ 880)  26.6 W49-W55 S08-S05
      9 1075  839 (S08W04 L= 93)   262    136        1 (1040/ 821-1117/ 870) -21.8 W01-W06 S09-S06
     10  915  702 (S17E06 L= 83)  1383    386        4 ( 804/ 582- 978/ 859)  68.7 E14-E02 S25-S07
     11  192  732 (S19E70 L= 20)   839    104        0 ( 164/ 690- 235/ 755)  15.5 E79-E61 S22-S17
     12  492  733 (S17E36 L= 53)   249     50        1 ( 472/ 709- 516/ 755)  45.0 E37-E33 S18-S15
     13 1230  680 (S18W14 L=103)   244     77        0 (1195/ 654-1257/ 713) -45.0 W11-W16 S19-S15
     14 1497  535 (S28W37 L=126)  1681    258        2 (1455/ 442-1543/ 656) -63.4 W33-W44 S35-S19
     15  349  541 (S31E56 L= 33)   607    109        0 ( 317/ 496- 382/ 572)  23.2 E64-E50 S35-S28
     16 1694  496 (S31W60 L=149)   563     79        1 (1654/ 485-1731/ 511)   6.3 W54-W65 S32-S30
  

    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