Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1468 1367 (N13W36 L=281)  1510    361        3 (1349/1303-1610/1395) -11.3 W26-W46 N08-N16
      2  190 1116 (N31E57 L=187)  1612    250        6 ( 158/1025- 249/1201)  68.2 E61-E51 N25-N36
      3  796 1058 (N15E12 L=232)  1852    477        4 ( 598/1026- 990/1111) -20.6 E24-W00 N09-N20
      4 1909 1066 (S18W63 L=308)   758    186        0 (1872/1041-1946/1087)  -9.5 W57-W69 S19-S17
      5 1819  903 (S24W49 L=294)   735    137        1 (1788/ 828-1838/ 973)  49.4 W45-W52 S28-S20
      6 1434  870 (S13W19 L=264)   398    150        0 (1407/ 798-1472/ 950)  45.0 W16-W23 S16-S09
      7  162  807 (N13E69 L=175)  1655    224        3 ( 119/ 749- 201/ 857) -40.6 E74-E66 N08-N15
      8  466  682 (N01E43 L=201)   489     70        1 ( 432/ 651- 498/ 708) -14.0 E45-E41 S01-N02
      9 1432  626 (S28W14 L=258)   268     45        0 (1407/ 601-1457/ 654)  11.3 W11-W16 S29-S26
     10  840  458 (S21E27 L=217)  1821    328        3 ( 731/ 309- 916/ 593)  79.7 E32-E25 S32-S10
     11  316  493 (S08E67 L=177)   577     56        0 ( 291/ 467- 340/ 521) -63.4 E69-E65 S10-S06
     12  688  338 (S26E43 L=201)   393    138        0 ( 667/ 297- 725/ 390)  69.4 E44-E40 S29-S21
  

    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