Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  509 1529 (N35E43 L=157)  1382    276        3 ( 408/1482- 627/1593)  -6.3 E58-E31 N31-N38
      2  711 1475 (N32E24 L=177)   261     59        0 ( 681/1453- 737/1508) -31.0 E26-E21 N30-N33
      3  762  686 (S18E17 L=184)   421    222        0 ( 666/ 661- 854/ 726) -13.0 E23-E10 S19-S15
      4  412  576 (S27E47 L=154)  2091    316        5 ( 325/ 426- 489/ 687)  30.3 E73-E37 S40-S19
      5  820  609 (S23E13 L=187)   338    176        0 ( 770/ 546- 890/ 651)  37.9 E17-E08 S27-S20
      6 1203  498 (S31W13 L=214)   265    104        0 (1159/ 476-1255/ 519)  23.2 W10-W17 S32-S29
      7 1590  427 (S37W51 L=251)   628     57        2 (1566/ 394-1611/ 459) -35.5 W46-W54 S39-S34
      8 1470  362 (S42W41 L=242)  1135    126        3 (1418/ 337-1536/ 395)  10.3 W36-W48 S44-S39
      9 1175  314 (S46W15 L=215)   330     49        1 (1155/ 293-1193/ 344)  63.4 W13-W16 S48-S43
     10 1269  304 (S47W24 L=225)   732    121        1 (1212/ 285-1332/ 328)   0.0 W18-W29 S49-S44
  

    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