Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  963 1595 (N38E11 L=299)   357     86        1 ( 928/1567-1006/1618) -18.4 E13-E07 N36-N39
      2 1111 1532 (N35W01 L=311)   278     66        0 (1082/1512-1142/1548) -18.4 W00-W03 N33-N36
      3  514 1455 (N24E41 L=269)   448    102        2 ( 482/1419- 545/1483) -33.7 E43-E37 N21-N25
      4 1266 1425 (N29W13 L=323)   249     59        1 (1244/1400-1296/1449)  26.6 W11-W15 N27-N30
      5 1444 1335 (N24W26 L=336)  1426    263        2 (1325/1248-1544/1408) -29.4 W17-W33 N19-N28
      6 1077 1241 (N15W02 L=311)   190     70        0 (1065/1200-1087/1277) -68.2 W00-W02 N12-N17
      7 1089 1145 (N10W03 L=313)   271     96        1 (1036/1124-1124/1181)  36.9 W00-W04 N08-N12
      8 1588 1091 (N10W37 L=346)   385     86        1 (1549/1056-1621/1128) -38.7 W33-W39 N08-N12
      9  147  966 (S10E69 L=241)   879    142        2 ( 132/ 891- 159/1035) -71.6 E71-E66 S14-S05
     10  553  750 (S19E29 L=280)   621    214        0 ( 473/ 709- 637/ 796)  30.3 E35-E23 S22-S15
     11 1837  731 (S11W66 L= 16)  1385    262        2 (1797/ 694-1883/ 768) -14.9 W59-W75 S13-S08
  

    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