Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1057 1598 (N34E12 L=183)  3161    702        5 ( 779/1467-1364/1747)  -5.7 E38-W12 N30-N39
      2  413 1602 (N21E59 L=136)   920    248        0 ( 392/1514- 431/1666) -45.0 E65-E54 N15-N25
      3  647 1442 (N15E33 L=161)   397     75        2 ( 620/1404- 682/1485)  90.0 E34-E31 N12-N18
      4  800 1354 (N13E20 L=174)   242     64        1 ( 778/1328- 815/1389)  76.0 E21-E19 N11-N15
      5  424 1270 (N01E44 L=151)   707    190        1 ( 361/1229- 520/1324)  35.0 E48-E37 S02-N05
      6  788  972 (S10E13 L=181)   427    133        0 ( 759/ 917- 840/1040) -40.6 E15-E08 S11-S05
      7 1885  990 (N15W64 L=258)   466     51        1 (1863/ 969-1905/1009)   0.0 W60-W66 N13-N15
      8 1016  928 (S08W01 L=196)   260     73        1 ( 973/ 899-1059/ 953)  53.1 W00-W03 S09-S05
      9 1694  667 (S08W54 L=248)   598    115        0 (1681/ 622-1710/ 720) -45.0 W51-W56 S10-S05
     10 1360  610 (S19W30 L=225)   234     49        0 (1334/ 595-1390/ 627)   0.0 W28-W32 S20-S18
     11 1038  554 (S30W12 L=206)   323     66        0 (1001/ 532-1069/ 573)  -9.5 W08-W14 S30-S28
  

    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